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Abstract 



m 
o 

I present a selection of CP-PACS and UKQCD data for the pseudo-Goldstone masses in 
Nf = 2 QCD with doubly degenerate quarks. At least the more chiral points should be 
consistent with Chiral Perturbation Theory for the latter to be useful in an extrapolation 
to physical masses. I find consistency with the chiral prediction but no striking evidence 
for chiral logs. Nonetheless, the consistency guarantees that the original estimate, by 
Gasser and Leutwyler, of the Nf = 2 QCD low-energy scale A3 was not entirely wrong. 
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Introduction 



Ever since numerical lattice QCD computations have been done, the spectrum of light mesons 



in 
o 

OO ' nas served as a benchmark problem. This is still true today, as actions respecting the global 
chiral symmetry among the light flavours are being tested and several groups have embarked 
on ambitious simulations of "full" QCD with two or possibly more dynamical flavours. These 
developments address one of the key issues in low-energy QCD. The fact that chiral symmetry 
is both spontaneously and explicitly broken generates pseudo-Goldstone bosons, i.e. particles 
that dominate (for small enough quark masses) the long-range behaviour of correlators between 
external currents and which are collectively called "pions" . 

The lattice is not the only framework to address the low-energy structure of QCD. In the 
old days PCAC relations were exploited to predict the dependence of low-energy observables 
on shifts in the quark masses and external momenta. The one best known is 

FlMl = (m u +m d ) |(0|gg|0)| + 0(m 2 ) (1) 

which connects the pion mass and decay constant to the product of the explicit and spontaneous 
symmetry breaking parameters. However, the Gell-Mann-Oakes-Renner relation does not 
give a prediction how F n and M v separately depend on the quark mass. And in the real world 
the latter may be shifted only by a discrete amount (e.g. by replacing d — > s,n — > K one gets 
a leading order prediction for the quark mass ratio m s /m with m = (m u +md)/2). 

Today, the first limitation is overcome, since the successor of PCAC, Chiral Perturbation 
Theory pQ, gives detailed predictions how either or F w individually depends on the quark 
mass m. And - in principle - the second problem is gone since one may vary m continuously 
in a lattice computation. Therefore, it seems natural to combine the two approaches to benefit 
from their respective advantages. However, for that aim quarks need to be taken sufficiently 
light, and this is a numerical challenge on the lattice. Below, an elementary test is presented 
whether this is already achieved in present state-of-the-art results for M 2 versus m, as published 
by the CP-PACS j2] and UKQCD j3] collaborations. The idea is to restrict the analysis to the 
"doubly degenerate" case with Nf = 2 0(a) improved Wilson quarks, i.e. to consider only the 
subset where both valence-quarks are exactly degenerate with the sea-quarks, as this avoids 
additional assumptions of the "partially quenched" framework. 
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Prediction by Chiral Perturbation Theory 



Here, I give a brief summary of the work by Gasser and Leutwyler (GL), with an application 
to lattice data in mind. 

GL have calculated the substitute for the GOR relation (pQ) to NLO in Chiral Perturbation 
Theory (XPT) with Nf =2 quarks pQ. To that order, two low-energy constants 1 from the LO 
Lagrangian appear (F = lim m ^ F n , B = — lim m ^ (0|gg|0)/F 2 ) together with the finite parts of 
the li in the NLO Lagrangian. Restricted to the degenerate case their result pQ takes the form 



Jif2 JUT2 

' logH^+M 2 ^)) (4) 



FV 16vr 2 to > 2 
M 2 = 2mB . (5) 

Note that (JHJ EJ) has the typical structure of a NLO prediction: A chiral logarithm summarizing 
the contribution from the pion loops appears together with a counterterm. The l\ are the 
renormalized GL coefficients, i.e. they are the descendents of the U which appear in the NLO 
Lagrangian and which are divergent quantities. As a result, the l\ depend on the (chiral) 
renormalization scale [i. In the dimensionally regularized theory one has pQ 

k = Z^ + ^AO) (6) 

w s -i ( i , iog(47r) + r(i) + i A (r . 

IIM = W-^log^) (8) 

and the /5-function coefficients 7^ are known. In the present context only 73 = — | , 74 = 2 pQ 
are relevant, and this suggests that one rewrites (0111) with the help of 

where the /i dependence in /J = /[(/i) is traded for an M dependence in li = li(M), to get pQ 

M 2 

Ml = M 2 (l - ^-^h + 0(M 4 )) (10) 
M 2 

F « = F ( 1 + T^ h + 0{M4) )- 

The NLO part is given in terms of the LO parameters F, B and the NLO coefficients h,U. 
While the former two are known quite accurately, for the U only their running in M 2 is known 
exactly and the phenomenological estimate of the integration constants has comparatively 
large error-bars. Gasser and Leutwyler give in their initial paper [Tj the estimate 

^ 3 (m phys ) = 2.9 ± 2.4 , U{m phjs ) = 4.3 ± 0.9 (12) 



1 Unlike F, the other LO "constant" B is scheme- and scale-dependent, as is m. In the chiral representation, 
only the product mB appears which is RG-invariant (cf. ©)■ In order to determine B and m separately, a 
lattice computation is needed; the result of the CP-PACS study along with © for the physical pion is 

m phys (MS,/i~2GeV) ~ 3.5±0.2MeV , B(MS, GeV) ~ 2.8 =F 0.15 GeV . (2) 
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for real world quark masses. It turns out that even this limited information is useful, since 
it determines the curvature of M n as a function of the quark mass. Close to the chiral limit, 
both li are positive and as a consequence M 2 does not rise strictly linear in m but turns right, 
while F n has a positive first derivative in m. More specifically, (|12j) translates into 

M^ phys ~139MeV, F 7r , phys ~92.4MeV <^=> M phys ~ 141 MeV , F~86.1MeV (13) 

which means that the physical pion is somewhat lighter than it would be, if the LO relation 
were exact, while is decay constant exceeds its value in the chiral limit by ~7%. 

A seemingly formal point which, in the end, proves convenient in analyzing the lattice data 
is the following. A naive look at (fTUJ) suggests that the typical structure of a NLO prediction is 
gone - rather than a M 4 and a M 4 log(M 2 ) contribution, only the polynomial part seems left. 
The point is that this impression is entirely misleading; the IR divergencies (which are genuine 
to QCD in the chiral limit) are not gone, they are just hidden in the Zj. The situation is, in 
fact, opposite - the M 4 part has been eliminated in favour of a pure M 4 log(M 2 ) contribution, 
and the last step is to make this apparent. The quark mass dependence of l 3 is given through 

r 2 ■ 

phys 

and together with (|5|). the relation takes its final form (see e.g. [3]) 



log(^) = Mj^) + M 1 ^) = ^Khys) + = U™) (14) 



M l = 2mB ~^-2 M^) + 0(m") (15) 



mB A 2 
h^F ° g ^2mB 

where A 4 ,A 3 represent universal low energy scales. The estimates ()12jl translate into 



F +^ lo g(^) + 0(m 3 ) (16) 



a n r~i +1.4 GeV . ^ „ „ T . +0.7 GeV 

A 3 = 0.6 GeV _ Q 4 ^ , A 4 = 1.2 GeV _ Q A ^ . (17) 

It is worth emphasizing that all four parameters F, M 2 = 2Bm, A 3 , A 4 do not depend on the 
QCD renormalization scheme pQ, i.e. they are the proper physical low-energy parameters at 
order 0(p 4 ). Furthermore the A; do not depend on the quark masses, and this means that the 
representation (|TK| Hfi]) is perfectly suited to analyze the lattice data even if they are gotten at 
quark masses larger than m = (m ll +m < i)/2 in the real world - as long as they are not beyond 
the regime of validity of the chiral expansion itself. 

The latter point is one of the key issues in the comparison we aim at. The chiral expansion 
is known to be asymptotic, and this means that increasing the order will enhance the accuracy 
near the chiral limit - at the price of worsening the prediction for heavier masses. What is the 
"critical scale" beyond which the chiral expansion "explodes" is, a priori, not known. From a 
formal point of view, one might think that (fTUl [TTj) indicate that the expansion is in Mj (4nF) 
and hence hope that it is good for pion masses up to 1 GeV. In this paper I will argue that 
watching the convergence pattern at a fixed quark mass gives a more reliable estimate what is 
the permissible range. This is facilitated since the NNLO expression for M 2 (with Nf =2) is 
known jHIISl- The result reads (for the presentation I follow 0]) 

9 ( mB , Ao m 2 B 2 f 17/, / A 2 » /fN \2 , , „ 
Ml = 2mB[ 1 — bgMr + -r^\—{ og(^ ) + k M \ + m 3 



A 2 A 2 A 2 A 2 

where A M is implicitly defined through 51 log(— f-) = 281og(-J-) + 321og(^) — 91og(-|) + 49 
and the mass-independent accounts for the remainder at 0(p 6 ), in particular the new 
counterterms. Phenomenological values for Am and hu will be mentioned below. 
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CP-PACS ((3=2.1 ) and UKQCD, N f =2, m sea =m vah pert, renormalization with c sw as used in simulations 
10 i , , 1 1 1 
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CP-PACS VWI, pert. ren.. 
CP-PACS AWI, pert. ren.. 
UKQCD VWI, pert. ren.. 
UKQCD AWI, pert. ren.. 
CP-PACS VWI, bare 
CP-PACS AWI, bare 
UKQCD VWI, bare 
UKQCD AWI, bare 
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CP-PACS ((3=2.1 ) and UKQCD, N f =2, m sea =m va |, pert, renormalization with c sw =1 throughout 
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Figure 1: CP-PACS and UKQCD data converted to the form M% vs. 2m with masses in units of 
Tq 1 . Top: Perturbative renormalization with csw as used in the simulations. Bottom: Same but 
c sw = 1 throughout. A quadratic fit, constrained to go through zero, is applied to the renormalized 
AWI data. Segments of the asymptotic slope in the chiral limit highlight the curvature. 

4 



Lattice Data 



We are now in a position to esteem the results by the CP-PACS and UKQCD collaborations 
for the quark mass dependence of the pion mass. We shall consider, out of these Nf =2 data, 
only the two-fold degenerate case where both valence quarks have the same mass and where 
they are, at the same time, exactly degenerate with the sea quarks so that the theory is unitary 

The CP-PACS collaboration has simulated various ((3, k) combinations with an RG im- 
proved gauge action and a mean- field improved clover quark action [2|. They use a grid of 
size 12 3 x 24, 16 3 x 32, 24 3 x 48, 24 3 x 48 at f3 = 1.8, 1.95, 2.1, 2.2, respectively, which leads to a 
lattice spacing, if determined through the p mass, between 0.215 fm and 0.087 fm and hence a 
spatial box size between 2.58 fm and 2.08 fm. With so much information at hand, one could, 
in principle, attempt a continuum extrapolation for M 2 versus the sum of the (degenerate) 
quark masses, 2m. However, non-perturbative renormalization might be necessary, and/or the 
lattice spacing might be too large at the lower (3 values. For this reason, I have decided to 
concentrate on the 13 = 2.1 data, since here discretization effects are not supposed to be too 
large, and good statistics is available. At that (3 value, even the lightest pion is unlikely to 
suffer from finite-size effects, since M 7T L>7, and the lattice spacing determined via the p mass 
is of order 0.1 fm and hence comparable with that in the UKQCD simulations. 

The UKQCD collaboration works on a 16 3 x 32 grid, using different actions: Wilson glue 
and non-perturbatively 0(a) improved clover quarks [3]. Another point in which they differ 
from CP-PACS is a tactical one: they try to relax (3 in pushing K sea up (i.e. the quark mass 
down) such that the lattice spacing, in units of ro, stays constant. Numerically, it is a~0.1fm 
for the data considered below (though the ensemble at (/3, n sea ) = (5.2, 0.1355) is not matched 
any more), and the hope is that the size of discretization effects would be approximately 
constant. This choice implies that the physical box size stays constant, too, and the bound 
M n L> 4.5 maintained makes one feel comfortable that finite size effects are small. 

The plan of this article is to ignore that the dynamical quark masses might be too heavy for 
the chiral prediction to be applicable and to ignore that in principle a continuum extrapolation 
is needed but instead to go ahead and simply compare the two datasets on a "as-is" basis with 
the LO/NLO/NNLO prediction from Chiral Perturbation Theory in the continuum. With the 
relevant low-energy constants on the chiral side given in physical units, Tq 1 must be so, too. 
In this article, this is done through the assumption that Tq 1 represents a universal low-energy 
scale, unaffected by unquenching effects; the numerical value used is r = 0.5fm [7j. 

On the lattice, there are two definitions of the "quark mass", one through the vector 
Ward-Takahashi identity (VWI mass), one through the axial identity (AWI mass). The bare 
masses (open symbols in Fig. 1) need not agree, while after renormalization (filled symbols 
in Fig. 1) they would (up to 0(a 2 ) effects) if the renormalization factors were computed non- 
perturbatively. For unquenched (Nf = 2) data the relevant non-perturbative factors are not 
yet available. For this reason, I have decided to renormalize both the CP-PACS and the 
UKQCD data at one- loop order (the details being given in the appendices). In the case of the 
CP-PACS data this basically repeats their calculation j2], albeit with two notable differences: 
First, the scale is set through the measured r , since this is the only possibility with the 
UKQCD data. Second, the "boost factor" uo is derived from the measured plaquette rather 
than the plaquette in the chiral limit, since in the UKQCD data there is no uniquely defined 
chirally extrapolated version. This perturbative calculation then lays the ground on which the 
CP-PACS and UKQCD data may be compared to each other and to the chiral prediction. 

Fig. 1 displays (M^tq) 2 versus the renormalized quark masses with filled symbols, open 
symbols indicate the bare data to visualize the shift. In the upper part the renormalization 
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CP-PACS (P=2.1) and UKQCD, N f =2, m sea =m va | 
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CP-PACS VWI, pert, ren.. 
CP-PACS AWI, pert. ren.. 

UKQCD VWI, pert. ren.. 

UKQCD AWI, pert. ren.. 

LO-XPT (B from phenomenology) 
LO-XPT, B -> B+1o 
LO-XPT, B->B-1cr 
NLO-XPT (B F A 3 from phenom.) 

N LO-XPT, A 3 ->A 3 -1o 
NLO-XPT, A 3 -> 
NNLO-XPT (A M from phen, ki 
NN LO-XPT 



-> A 3 +1a 



NNLO-XPT, k 




Figure 2: LO/NLO/NNLO chiral predictions with phenomenological values for B, F, A3, Am, &m 
(±1(7 variation included at each order, cf. text for details) compared to the renormalized data in the 
version with csw = 1- Note that the lines represent parameter-free predictions, not fits. 



was performed with the csw values as they were used in the simulations (suggested by a mean- 
field analysis [2. or the Alpha study 0IB!)- In the lower part the renormalization factors were 
computed with csw = l throughout, which is a consistent choice at one- loop order. Obviously, 
the overall consistency of the data is much better with this latter choice, as is particularly 
obvious from the quadratic fits (constrained to go through zero) to the AWI data. Henceforth, 
we shall stick to the latter choice (csw = l), but it is useful to keep in mind that the difference 
to the upper part is a clear indication of the size of inherent perturbative uncertainties. Finally, 
it is worth mentioning that all error-bars in this article represent only statistical errors. 

We now turn to the physics content. To convey a feeling for the scales, I mention that at 
2mro = 1 the sum of the valence quark masses is of order 400 MeV, i.e. about four times as 
much as in a physical K, and the corresponding "pion" weighs about 1 GeV. The physical 
kaon weighs 496 MeV, i.e. (MkTq) 2 ~ 1.58. If the pion would satisfy (M n ro) 2 — 1.58 this would 
mean that its u- and (i-quarks are about half as heavy as the s-quark in the real world. The 
lightest pion in the CP-PACS and UKQCD simulations have (M n r ) 2 values 1.75 and 2.20, 
respectively, and from this we conclude that their lightest u- and d- quarks (in the unitary 
theory) have about 55% and 70% of the mass of the physical strange quark. 

With such heavy "light" quarks it is a priori not clear whether XPT is of any use to 
extrapolate them to the physical u- and c/-quark masses. In an attempt to shed a light on 
this issue, Fig. 2 shows the renormalized data (as in the bottom part of Fig. 1) together 
with the predictions from XPT at tree/one- loop/two-loop level (LO/NLO/NNLO). The low- 
energy constants are taken from phenomenology, i.e. these curves represent parameter-free 
predictions. At LO the chiral prediction is a straight line with the slope parameter B taken 
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relative shift at LO/NLO/NNLO level of XPT w.r.t. previous order 
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Figure 3: Relative shifts in M% versus 2mro in XPT - due to the uncertainty in B at LO (light), 
due to NLO contribution with B fixed at central value, A3 varied within la hound (intermediate), 
and due to NNLO contribution with B and A3 fixed at central values, Hm £ {0, ±2} (dark). 

from (J2J); the ±lcr bounds are indicated by dotted lines. At NLO the prediction is curved, and 
the numerical values of the additional parameters are taken from (fT3*|) and (fTTj) . A3 is varied 
within its phenomenological ±ler bound (full versus dotted lines - lowering A3 to 200 MeV 
yields the upper, increasing it to 2 GeV yields the lower dotted line), with B, F fixed at their 
central values. At NNLO the parameters Am, kht in ()18|) are determined as follows. With 
Ai = 0.11±J;JJ GeV, A 2 = 1.2±0.06 GeV [9 and jfjj the relation beneath eqn. flU)) gives 

A M = 0.60±0.03GeV (19) 

where errors have been added in quadrature. This means that over the range considered 
the uncertainty in Am is neglibigle compared to that coming from ku- Phenomenological 
arguments indicate \ku\ ~ 2 ^0], and the sum- rule estimates for the NNLO counterterms 
given in [H] may be converted into a more accurate estimate. For the purpose of the present 
article it is sufficient to use k M = ± 2, and the associate curves (holding B,F,A 3 ,A M fixed 
at their central values) are included in Fig. 2. 

This information is sufficient to asses the chiral convergence behaviour in M% versus 2m. To 
that aim we compare the LO/NLO/NNLO predictions for a given quark mass. Fig. 3 shows the 
relative shift in M% when one more order is included or the new low-energy constants are varied 
within reasonable bounds. At first sight, the uncertainties at higher orders due to the error 
bars of the counterterms seem large, but one should keep in mind that the associate shifts are 
100% correlated over the whole range. For instance, if the CP-PACS VWI point at 2mr ~0.2 
sits on the —la curve (upper dotted NLO line in Fig. 2), then the one at 2mro~0.4 should 
too, if only A 3 needs to be adapted. This means that precise lattice data are ideally suited 
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Figure 4: Same as Fig. 2, after dividing by 2mrQ. Still, lines represent parameter-free chiral predic- 
tions at LO/NLO/NNLO (from light to dark colour), not fits. 



to reduce the error on ll(fi) (or A3). Fig. 3 shows that the crossover where the uncertainty 
due to the NLO exceeds that due to the LO contribution is at 2mro = 0.45/0.02/0.10 when 
A 3 is taken at its central/ + la/- la value. By averaging with weights 2/1/1 one arrives at the 
estimate that in the case of M x the chiral expansion is sufficiently well-behaved that the NLO 
functional form is useful for quark masses up to 

2mr < 0.25 (M n r ) 2 < 2 M^< 560MeV, (20) 

and maybe more. This, if correct, means that current state-of-the-art simulations make contact 
with the regime where NLO-XPT holds, but so far a non-trivial "lever arm" which is needed 
to make model independent predictions in the deeply chiral regime, is likely missing. 

While it is clear that future simulation data will allow to test the prediction (j2*Uj) and 
provide, if it is correct, the lever arm needed, one might, already at this time, go ahead and 
try what comes out if one assumes that the estimate (|2T)|) is too pessimistic and hence uses the 
NLO chiral ansatz to fit the data over an extended range. This is what we shall do below, but 
it is clear that this attempt is rather speculative and results should be taken with care. 

Since everything below is about the deviation from a linear relationship, it is useful to make 
the curvature optically visible. This is conveniently done by dividing out a factor 2mro- The 
result is displayed in Fig. 4 where the predictions from XPT at LO/NLO/NNLO are included 
for completeness. In this representation, the genuine feature of the NLO curve is that it lies 
below the LO constant for light pions, but above if the (LO-) pion mass is larger than A3. 

To fit the data, one replaces 2mB by M%, since in this form the data are plotted against 
something which is directly measured and the change is of yet higher order in the chiral 
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Figure 5: Fits of the data with the AWI quark mass to the NLO functional form (together with the 
allowed band from phenomenological estimate of A3). In spite of the bound i|20j) being ignored, the 
data ht well to the suggested form (even if one insists on the phenomenological value of F n in the 
NLO ansatz i)151 \21\) ). but with the current error-bars the LO functional form is still appropriate. 
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8 05 f0 234") 3 21 73/2 
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10 fO 234") 3 95 04/2 
9.97 (0.234) 3.93 0.04/1 
10.1 0.231 3.94 0.04/1 


UKQCD-LO 


6.67 1.39/3 


6.90 1.51/3 



Table 1: Coefficients in the fits of the NLO functional form ^71} to the degenerate CP-PACS and 
UKQCD data with VWI (left) or AWI (right) definition of the quark masses, after 1-loop renor- 
malization as detailed in the appendix. Phenomenological values for comparison; constrained values 
in brackets. In one case, there is such a shallow minimum that the fitting algorithm fails to determine 
it. In general, x 2 is underestimated, since the fact that errors are correlated has been ignored. Since 
the main uncertainty is systematic, I refrain from quoting statistical errors. 

counting. In other words, the statement is that one should use the representation 

^^{l-^Llog(S)} (21) 

and adjust the 3 dimensionless parameters B = Br , F = F n r , A 3 = A 3 tq. 

Attempts to fit the data with either the VWI or the AWI definition of the quark mass to 
formula (}2*T]) are summarized in Tab. 1 and - for the AWI mass - illustrated in Fig. 5. Since 
the parameter F n = F n ro = 92.4 MeV 0.5 fm ~ 0.234 is known, it is held fixed at its physical 
value in (a, b) and fitted only in case (c). Fit (b) differs from (a) in that the data point with 
the heaviest quark mass has been omitted. From the overall spread it is clear that there are 
considerable systematic uncertainties and this is why I refrain from quoting statistical errors. 

It is clear that with the quality of the present data one cannot claim direct evidence for 
chiral logs. Nonetheless, it is worth noticing that fits (a) and (b) indicate that the data 
are consistent with the logarithmic form (}2*T]) suggested by NLO chiral perturbation theory. 
Moreover, that the fit with F v held fixed at its physical value does not dramatically change if 
the heaviest data point is omitted (a— >b), suggests that the estimate (f2U|) for the permissible 
range at NLO level is indeed reasonable. What one can learn from Table 1 is that in QCD 
with Nf = 2 the low-energy parameter A 3 does not differ by orders of magnitude from the 
original estimate by Gasser and Leutwyler pQ; all A 3 values in Tab. 1 lie between the central 
and the +lcr level in (|T2*|) . This is important, since it excludes - at least for Nf = 2 - an 
alternative scenario of the chiral symmetry breaking in which B would not be the adequate 
order parameter jTTj (for an introduction to this topic see j3]). 

Discussion 

The aim of the present note has been to compare the CP-PACS and the UKQCD data for M% 
as a function of the quark mass to each other and to the prediction from Chiral Perturbation 
Theory (XPT). We have seen that after renormalizing with 1-loop lattice perturbation theory 
the two sets are reasonably consistent (maybe with the exception of the CP-PACS data with 
the VWI definition of the quark mass). After dividing out a factor 2mr to "zoom in" on the 
deviation from a linear behaviour, both sets show very little curvature. Postponing the issue 
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whether or not it is permissible to compare non-continuum-extrapolated lattice data to the 
prediction from XPT in the continuum we found that the current size of the error-bars assures 
consistency with the logarithmic form (J21|) . albeit one cannot claim evidence for chiral logs. 

The situation might change as more precise data are being released. In fact, the JLQCD 
collaboration claims that their data with Nf = 2 non-perturbatively 0(a) improved clover 
fermions and standard (Wilson) glue are inconsistent with NLO chiral perturbation theory, if 
they keep the parameter F n fixed at its physical value ^2j . 

At this point it is mandatory to think about possible reasons why the data might not 
coincide with the chiral prediction. An incomplete list is the following one: 

1. They need not, if the pions are too heavy. Most of the data have been obtained in a 
regime of quark masses where there is no guarantee that XPT works. Phenomenological 
experience tells us that this expansion works for mesons as heavy as the physical K, i.e. 
for (M n r ) 2 < 1.58. I have argued that this condition might be relaxed to (J2"U|) . Even if 
this is true, at best 1 point from either set survives, and hence there is, strictly speaking, 
no room for trying the NLO-XPT ansatz which has 3 parameters. 

2. Scaling violations, in particular implications of the broken chiral symmetry might be so 
severe that the chiral logs - even if they exist in the mass range considered - might be 
buried. To check one would have to perform a continuum extrapolation or attempt a 
dynamical simulation with fermions which obey the Ginsparg- Wilson relation. 

3. For unknown reasons (e.g. an algorithmic flaw), the data might represent a partially 
quenched rather than the fully unquenched situation. Obviously, this is a rather remote 
possibility. The point is that the NLO prediction M 2 /M 2 = 1 + constM 2 log(M 2 /A^) 
has no genuine M 2 term. On the other hand, in the partially quenched case a true M 2 
contribution exists [T4*| I15j: hence even perfectly linear behaviour in a plot analogous to 
Fig. 4 or 5 does not necessarily imply unreasonable values of the low-energy constants. 

Finally, I would like to highlight 4 points: 

(i) So far, we have ignored systematic uncertainties. Going back to Fig. 1 one sees that they 
are far from being negligible. What we see exemplifies the standard wisdom that perturbative 
renormalization factors which turn out to deviate from 1 by, say, 10% call for either repeating 
the exercise with non-perturbatively determined Z-factors or at much smaller lattice spacing. 
In the present context this means that the figures in Tab. 1 should not be taken as lattice de- 
terminations of the low-energy constants F,B,A^. Nonetheless, even this preliminary attempt 
shows that A 3 cannot differ by orders of magnitude from the original estimate by Gasser and 
Leutwyler [I]. Moreover, as the data get more precise - and under the proviso that consis- 
tency is maintained - lattice determinations of the scales A 3 and A4 in conjunction with chiral 
perturbation theory will yield precise predictions of intrinsically Minkowskian quantities, e.g. 
the I = 0, 2 pi-pi scattering lengths (see 0] for the connection). 

(ii) For a thorough comparison with chiral perturbation theory, one needs to get control 
over the lattice artefacts. A safe way to do this is to first extrapolate all data to the continuum. 
Then one can determine the permissible mass range and extrapolate in a second step to the 
chiral limit. If it turns out that this order cannot be sustained in practice, the chiral framework 
may be extended to account for the main discretization effects. For the case of unimproved 
fermions this has been done and it is clear that this approach could be generalized to 
improved actions. The only disadvantage is the added number of counterterms that need to 
be fixed from the data. 
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(Hi) In the present note the scale is set through Tq 1 throughout. In the CP-PACS studies, 
f has been found to strongly depend on m VWI or Ej. If r itself depends on the 
sea-quark mass, normalizing all masses through the measured f may not be the best way to 
compare to XPT; maybe one should use f 0i o (the chirally extrapolated version) instead. 

(iv) Strictly speaking, the "phenomenological" value of l T 3 (or A3) that we used in comparing 
the lattice data to the prediction from XPT is not quite adequate, since the phenomenological 
determination is with m s fixed at its physical value, while the lattice studies are Nf = 2 
simulations, in other words here m s is sent to infinity. The connection to the Nf = 3 GL 
coefficients L\(^j) is given by [IB] 

11 = 8(2L; - LI) + 4(2L' 8 - Li) - ^ ( log (^=^1) + l) (22) 

[the limit on the r.h.s. refers to a situation with m s held fixed at its physical value], and 
it is tempting to use it as the starting point of a little gedanken experiment: Assume the 
s-quark mass would be such that even when it is doubled the chiral expansion would not 
break down. Eqn. (J22|) tells us that under these circumstances enhancing m s by a factor 2 
(and hence roughly doubling M^, too) would lower l T 3 by 0.0001. Since (fT2|l translates into 
/!j(yU ~ M p ) = 0.0008 ± 0.0038 such a shift would be negligible compared to the error-bar, 
and it seems therefore reasonable to assume that the effect of an added third flavour in the 
simulations, fixed at the physical s-quark mass (which would fully justify the comparison with 
Nf =2 XPT), would be small compared to the inherent theoretical uncertainties. 

Formula ()22|1 is interesting in yet another respect, since it tells us that lattice studies which 
pin down M% versus 2m (and hence A 3 or l T 3 ) have a say in another issue. In the 3-flavour theory 
there is the famous "Kaplan Manohar ambiguity", i.e. the chiral Lagrangian stays invariant 
under a simultaneous transformation of the quark masses 

m u — > m n + Amdin s , rnj — > rnj + Am s m u , m s — > m s + \m u md, (23) 

and an appropriate modification of L r 6 ,L\,L T s . This would, in principle, allow to tune m u = 
and hence provide a simple and elegant solution to the strong CP problem JH]- in XPT 
terminology it is the low-energy combination 2L\((i) — L r 5 (fi) that decides whether this is a 
viable option [20] • In the past, XPT has been augmented by theoretical assumptions or model 
calculations to exclude m u = 0. More recently it has been proposed to determine the L\ from 
lattice simulations [21] and important steps in this program have been taken [22]- Formula 
([22]) tells us that a lattice determination of l\, augmented by knowledge about L\ and L T 6 , 
helps to constrain m u = 0. 

The present analysis certainly emphasizes the need to compute renormalization factors 
non-perturbatively and to perform a continuum extrapolation with dynamical data. 
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Appendix A: Renormalization of the VWI quark mass 

Setup with m VWI a = log(l + §(± - £)) ~ §(± - and u = P 1 / 4 = (i^)) 1 /* [23]: 

^SpG") = ^m(A*o)(l + 6 m am)m (24) 



= f m M(l + 6 m — )- (25) 
From [21], naive versus tadpole-improved, with g 2 = 6/(3 and g 2 = £feg(2 GeV): 

Z m (H = l + |!(^-Ilog(aV)) (26) 

Z m (H = i + £(^_Ilog(aV)) (27) 
From [25] (key to orig. literature: their [21-28]), for generic actions: 

\ + d g + d f N f + n ~ 2i ! //3 log(*V) (28) 



Wilson/Clover: d g = -0.4682, d f = 0.0314917 (for c sw = 1) and P = 1 - 1/3 • g 1 

1 1 11 - 2JV//3, , n o, 



- 0.4682 + 0.0314917iV) + J/ log(aV) (29) 



1 P 11 - 2N f /3, , 9 

0.1349 + 0.0314917iV) + /' log(aV) (30) 



z m = 12.953 + 7.738c sw - l-380c| w (31) 
z - z - vr 2 - { 13 -° (Csw " 2) (32) 

Iwasaki/Clover: d 9 = 0.1000, d/ = 0.0314917 (for c sw = 1) and P = 1 - 0.1402# 2 , 

P = 1 - 0.2689# 2 , thus 1 = P + 0.1402c/ 2 = P + 0.2689# 2 = 3.648P - 8 • 0.331P - 0.2006# 2 

1 1 11- 2M/3, , - 9N 
-5-TT = + 0.1000 + 0.0314917iV) + —f 1 log a 2 /i 2 33 

^Ms(^) 9 167r 

^ + 0.2402 + 0.0314917^ + - "7 ' " log(a'V) (34) 



1 p 11 O AT /O 

To — / — r = + 0.2402 + 0.0314917iV) + — -J' < 2 -> 

^Mg(^) 9 2 16tt^ 

IP 11 - 2N f /3 , ., 



+ 0.3689 + 0.0314917iV) + J/ log(aV) (35) 



1 _ 3.648P - 2.648P „ M( ,,, m ,- V LI -2^/3, 



0.1006 + 0.0314917iV) + " log(aV) (36) 

SKEW » 167r 

2 m = 4.858 + 5.301c sw - 1.267c| w (37) 

~z m = z m -0A206n 2 = { *™ ( ' SW = L ^ (38) 

4.74 (c sw = 1) 



In one finds: 



Wilson/Clover: b m = -1/2 - 0.09623# 2 + 0(# 4 ) (39) 

b m = b m (l-l/12-g 2 ) ~ -1/2 - 0.05456^ (40) 

Iwasaki/Clover : 6 m = -1/2 - 0.0509# 2 + 0(g 4 ) (41) 

b m = 6 m (l - 0.03505*? 2 ) ~ -1/2 - 0.0334<? 2 (42) 
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Appendix B: Renormalization of the AWI quark mass 

Setup with m AWI ' imp a = \(d^Af™v{x)0 & {U)) / (P & {x)0{$)) and u = P 1 / 4 = (|(M D ))V 4 p 



AWI^ \ ^ 1 + bAdm AWIjmp /ao\ 

(V) = ; — r ; m ' p 43 



m^(fi) = J A 1 + bAam/U ° m AWI ' imp 
Zp(fia) l + bpam/iiQ 


(AA\ 


J-IUIII IZi^l . lldlV" VCIoLlo LdLljJUlt 


J - iiii]jiuvcu. ; wiLii ty — u / diiu. y — ^ MS ^ * /• 




* ~ 1 + C 


^, = l + f + -MaV)) 


(45) 




MH = i + f (f + -MaV)) 


(46) 


From j27], for generic actions: 






Wilson/Clover: 






za = 


= -15.797 - 0.248c sw + 2.251c^ w 


(47) 


za = 


, / 2.581 (c sw -2) 
1 -3.924 (c sw = 1) 


(48) 


zp = 


: -22.596 + 2.249c sw - 2.036c| w 


(49) 


zp = 


, / -16.372 (c sw ^ 2) 
P+ \ -12.513 (c sw = l) 


(50) 


Iwasaki/Clover: 






z A = - 


8.192 - 0.125c sw + l-610c| w 


(51) 


, n ,on. 2 f -0-746 c sw = 1-47 
[ -2.556 (c S w = 1) 


\ oz ) 


zp = - 


10.673 + 1.601c sw - l-281c| w 


(53) 


2 / -6.936 (c S w = 1.47) 
z P = * P + 0.4206* = | _ 6 2Q2 (cgw = 1) 


(54) 


In one finds: 






Wilson/Clover : 


b A = 1 + 0.15219(5)/ + 0(g 4 ) 


(55) 




b A = b A (l-l/12-g 2 ) ~ 1 + 0.06886(5)/ 


(561 




6 P = 1 + 0.15312(3)/ + 0{g 4 ) 


(57) 




bp = 6 F (1 - 1/12- <f) ~ 1 + 0.06979(3)/ 


(58) 


Iwasaki/Clover : 


6 A = 1 + 0.0733(5)/ + 0(g A ) 


(59) 




b A = 6^(1 - 0.03505/) ~ 1 + 0.0383(5)/ 


(60) 




bp = 1 + 0.0744(12)/ + 0(/) 


(61) 




bp = b P {l - 0.03505/) ~ 1 + 0.0394(12)/ 


(62) 


In all cases the relationship between Zx and zx as well as bx and &x {X £ A P}) 


is mine. 
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Appendix C: Implementation with tadpole resummation 



The renormalization with uq defined via the plaquette is traced in Tab. 2 for the UKQCD 
data. Alternatively, one could define uo as 1/(8k c ), but this would bring different perturbative 
coefficients than those listed in App. A/B. Notice the effect of the tadpole resummation, i.e. the 
difference of the coupling computed via (|2*9*j) an d g^ via (J3U|) . The statistical uncertainty 
for the VWI quark mass is larger than that for the AWI definition due to a limited accuracy 
of k c , which is defined in a partially quenched sense 



UKQCD 



(5.20,0.1355) (5.20,0.1350) (5.26,0.1345) (5.29,0.1340) 



fo = r /a 
(M,r ) 2 



13 
13 



0.294(4) 
5.041(40) 
2.20(09) 



0.405(4) 
4.754(40) 
3.71(14) 



0.509(2) 
4.708(52) 
5.74(17) 



0.577(2) 
4.813(45) 
7.71(20) 



K c 

m 



VWI 



m 



AWI,imp 



1(1 _ ±) 

2 V K Kc' 



p |2H1 

(a2GeV) 2 = (5.06773/f ) 2 
^(2 GeV) m 
&(2 GeV) m 

bp fU) 



0.13645(3) 
0.0257(08) 

0.0231(3) 
0.536294(9) 

1.011(16) 
2.1640(45) 

2.541(6) 
-0.6387(3) 

1.1750(4) 

1.1773(4) 



0.13663(5) 
0.0442(13) 

0.0462(3) 
0.533676(9) 

1.136(19) 
2.1309(46) 

2.510(6) 
-0.6369(4) 

1.1728(4) 

1.1752(4) 



0.13709(3) 
0.0702(08) 

0.0742(3) 
0.539732(9) 

1.159(26) 
2.0813(58) 

2.437(8) 
-0.6330(4) 

1.1678(5) 

1.1701(6) 



0.13730(3) 
0.0897(08) 

0.0952(3) 
0.542410(9) 

1.109(21) 
2.0714(49) 

2.423(7) 
-0.6322(4) 

1.1669(5) 

1.1691(5) 



csw 

ZA 
Zp 

Z m (a2GeV) 
Z A 

Z P (a2GeV) 



dsn 

dm ehj 



2m^ 1 (2GeV)r W 
2m^ I (2GeV)r Q 



2.0171 1.9497 1.9192 

13.0769 12.9243 12.8512 

2.73099 2.14587 1.88782 

-16.4738 -16.081 -15.9094 

1.2799(17) 1.2690(18) 1.2569(22) 1.2566(18) 

1.0586(1) 1.0579(1) 1.0442(1) 1.0386(1) 

0.6472(1) 0.6590(2) 0.6781(3) 0.6808(3) 



0.380(16) 
0.381(09) 



0.603(25) 
0.705(13) 



0.919(23) 
1.076(19) 



1.181(24) 
1.398(21) 



csw 



z A , zp 



iny a2GeV) 



Z A 

Z P (a2GeV) 



SHI EH 

(Hi 



2mg 1 (2GeV)r JST 
2m^ I (2GeV)r Q 



1.2019(15) 

0.9158(2) 

0.7322(4) 



9.4414 , -3 
1.192(16) 
0.9168(2) 
0.7429(4) 



1 

9244 , -12.5134 
1.1852(20) 
0.9192(3) 
0.7516(5) 



1.1869(17) 

0.9197(2) 

0.7503(5) 



0.357(15) 
0.291(07) 



0.567(23) 
0.542(10) 



0.867(21) 
0.854(16) 



1.115(23) 
1.123(17) 



Table 2: Step-by-step renormalization of UKQCD quark masses with uq defined via the plaquette. 
Error bars include statistical errors only, naive error propagation throughout. Note the disparity of 
9mq defined via and that via if^D)) . It makes a difference whether csw as used in the simulations 



is plugged in or csw 



1 which is a consistent choice at the order we are interested in. 



For the CP-PACS data, both the plaquette P and the rectangle R are published [2], and 
this gives, in principle, the option to define the tadpole resummation (besides the usual option 
l/(8/c c )) via (j33j) or ([36)1 . where the latter choice reflects the specific combination used in the 
action. All these options would, however, imply different perturbative coefficients than those 
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listed in App. A/B, i.e. we restrict ourselves to the choice (JMj) along with uq = P 1 / 4 . It is 
interesting to see that the CP-PACS RG improved action achieves agreement of with the 
"standard" gLg (via flUD). 

In the argument of the logarithm that converts to the MS scheme, the lattice spacing is 
multiplied with a physical scale, here 2 GeV. This means that a must be assigned a physical 
value, as well. To compare like with like this is done via the measured ro in both sets, assuming 
tq = 0.5 fm in physical units. As a consequence, our Z factors for the CP-PACS data depend 
slightly on the quark mass, even though we work in a mass-independent scheme. 



CP-PACS 



(2.10,0.1382) (2.10,0.1374) (2.10,0.1367) (2.10,0.1357) 



f = r /a 
(M^r ) 2 



2J 

121 



0.29459(85) 
4.485(12) 
1.746(19) 



0.42401(46) 
4.236(14) 
3.226(28) 



0.51671(67) 
4.072(15) 
4.427(44) 



0.63010(61) 
3.843(16) 
5.864(60) 



rh 
m 



vwi _ . 

AWI,imp 



1(1 -±) 



P 

R 

(a 2 GeV) 2 = 
^GeV) 

?Ms( 2GeV ) 
~b A 



2J 

21 
21 
21 



(5.06773/fo) 2 
(ESI) 
(El 



0.02041(34) 
0.02613(18) 
0.6010819(84) 
0.366883(13) 
0.880508(3) 
1.277(07) 
1.8941(12) 
1.8920(12) 
1.7383(10) 
2.4705(20) 
-0.56315(4) 
1.07237(4) 
1.07445(5) 



0.1389; 
0.04147(34) 
0.05267(22) 
0.6000552(67) 
0.365297(10) 
0.880132(2) 
1.431(09) 
1.8694(14) 
1.8686(14) 
1.7191(12) 
2.4276(24) 
-0.56236(5) 
1.07147(5) 
1.07353(6) 



54(13) 
0.06011(34) 
0.07564(38) 
0.5992023(76) 
0.363979(12) 
0.879819(3) 
1.549(11) 
1.8526(15) 
1.8529(15) 
1.7063(13) 
2.3986(26) 
-0.56184(5) 
1.07087(6) 
1.07291(6) 



0.08706(34) 

0.10748(51) 
0.5980283(76) 
0.362139(12) 

0.879388(3) 
1.739(14) 
1.8286(17) 
1.8303(17) 
1.6877(14) 
2.3574(28) 

-0.56108(6) 
1.07001(7) 
1.07202(7) 



csw 

z mi z Ai Z P 

Z m (a2GeV) 
Z A 

Z P (a2GeV) 



2m^ i (2GeV)r 



(HID 
(Hi) 

dm) 



1.47 

5.76145 , -0.745545 , -6.93649 

1.08033(31) 1.07393(37) 1.06960(40) 1.06338(44) 

0.98809(1) 0.98824(1) 0.98834(1) 0.98848(1) 

0.90090(19) 0.90753(24) 0.91201(27) 0.91846(31) 



is 

2 mff I (2GeV)r 



(123) 



0.2217(44) 
0.2571(27) 



0.4174(50) 
0.4859(42) 



0.5723(57) 
0.6675(65) 



0.7642(66) 
0.8889(89) 



csw 

Zm-, %Ai Z P 

Z m (a2GeV) 
Z A 

Z P (a2GeV) 



2m^ i (2GeV)r 
2mM I (2GeV)r 



(J2ZD 
(gD 



1.06403(30) 
0.95917(3) 
0.91263(20) 



1 

4.74084 , -2.555^ 
1.05783(36) 
0.95968(3) 
0.91912(25) 



14 , -6.20184 
1.05363(39) 
0.96001(3) 
0.92351(28) 



1.04761(43) 
0.96050(4) 
0.92981(32) 



(123) 



0.2183(43) 
0.2463(26) 



0.4111(49) 
0.4659(40) 



0.5637(56) 
0.6403(63) 



0.7529(65) 
0.8532(86) 



Table 3: Step-by-step renormalization of CP-PACS quark masses with uq defined via the plaquette. 
Error bars include statistical errors only, naive error propagation throughout. Note the similarity 
of defined via ip£?j) and that via 1341). There is little difference whether csw as used in the 
simulations is plugged in or csw = 1 which is a consistent choice at the order we are interested in. 
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